Associated production of prompt photons and heavy 
quarks in off-sheh gluon-gluon fusion 

S.P. Baranov", A.V. Lipatov^ N.P. Zotov^ 
May 17, 2008 



P.N. Lebedev Physics Institute, 
119991 Moscow, Russia 

^ D. V. Skobeltsyn Institute of Nuclear Physics, 
M. V. Lomonosov Moscow State University, 
119991 Moscow, Russia 



Abstract 

In the framework of the /cT^-factorization approach, we study the production of prompt 
photons associated with heavy (charm and beauty) quarks in hadron-hadron coUisions at high 
energies. Our consideration is based on the amphtude for the production of a single photon 
associated with a quark pair in the fusion of two off-shell gluons. The total and differential 
cross sections are presented and the conservative error analysis is performed. Two sets 
of unintegrated gluon distributions in the proton have been used in numerical calculation: 
the one obtained from Ciafaloni-Catani-Fiorani-Marchesini evolution equation and the other 
from Kimber-Martin-Ryskin prescription. The theoretical results are compared with recent 
experimental data taken by the CDF collaboration at the Fermilab Tevatron. Our analysis 
extends to specific angular correlations between the produced prompt photons and muons 
originating from semileptonic decays of the final charmed or beauty quarks. We point out 
the importance of such observables, which can serve as a crucial test for the unintegrated 
gluon densities in a proton. Finally, we extrapolate the theoretical predictions to the CERN 
LHC energies. 

This paper is dedicated to the memory of P.F. Ermolov, who died on May 14, 2008. 
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1 Introduction 



The production of prompt (or direct) photons in hadron-hadron coUisions at the Teva- 
tron is a subject of intense studies on both the theoretical and experimental sides [1-9]. 
Usually the photons are called "prompt" if they are coupled to the interacting quarks. The 
theoretical and experimental investigations of such processes have provided a direct probe 
of the hard subprocess dynamics, since the produced photons are largely insensitive to the 
effects of final-state hadronization. At the leading order of perturbative Quantum Chromo- 
dynamics (pQCD), prompt photons can be produced via quark-gluon Compton scattering 
or quark-antiquark annihilation and so, the cross sections of these processes are strongly 
sensitive to the parton (quark and gluon) content of a protoE0. The perturbative QCD 
calculations [11] in the next-to-leading order (NLO) approximation agree with the recent 
high-pr measurements [6] within uncertainties (see also discussions in [12-16]). However, 
there are still open questions. It was found [1-4] that the shape of the measured cross sec- 
tion as a function of photon transverse energy Et is poorly described by the NLO pQCD 
calculations: the observed Et distribution is steeper than the predictions. This shape dif- 
ference leads to a significant disagreement in the ratio of the cross sections calculated at 
different center-of-mass energies \/s = 630 GeV and y/s = 1800 GeV as a function of scaling 
variable xt = 2Ej'/y/s. It was demonstrated [2, 3] that the disagreement in the xt ratio 
is difficult to explain with the conventional theoretical uncertainties coming from the scale 
dependence and different parametrizations of the parton distributions. In the NLO QCD 
approximation, the observed discrepancy can be reduced [13, 17] by attributting some addi- 
tional intrinsic transverse momentum to the incoming partons, which is usually assumed 
to have a Gaussian-like distribution. The average value of this increases from (kx) ~ 1 
GeV to more than (kx) ~ 3 GeV [13, 16] as the ^/s increases from UA6 to Tevatron energie^. 

From our point of view, a more adequate solution was found [22-25] in the framework of 
the /cy-factorization approach [26, 27]. This approach is based on the familiar Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) [28] or Ciafaloni-Catani-Fiorani-Marchesini (CCFM) [29] gluon evo- 
lution equations and takes into account the large logarithmic terms proportional to Inl/x. 
It is believed that such terms give a significant contribution to the heavy quark produc- 
tion cross section at the conditions of modern colliders. This contrasts with the usual 
Dokshitzer-Gribov-Lipatov-Altarelli-Parizi (DGLAP) [30] strategy where only the large log- 
arithmic terms proportional to In//^ are taken into account. The basic dynamical quantity 
of the fc^-factorization approach is the unintegrated (i.e., kT-dependent) gluon distribution 
A{x, h"^, fi"^) which determines the probability to find a gluon carrying the longitudinal mo- 
mentum fraction x and the transverse momentum k^" at the probing scale /i^. Similar to 
DGLAP, to calculate the cross sections of any physical process the unintegrated gluon den- 
sity A{x, k^, fj,"^) has to be convoluted [26, 27] with the relevant partonic cross section which 
has to be taken off mass shell (k^-dependent). Since the gluons in the initial state are not on- 
shell and are characterized by virtual masses (proportional to their transverse momentum), 
it also assumes a modification of their polarization density matrix [26, 27]. In particular, the 

"'^Also, the observed photon may arise from the so called fragmentation processes [10]. This contribution 
will be discussed below in Section 2. 

^The importance of including the gluon emission through the resummation formalism was recognized and 
only recently this approach has been developed [17-21] for inclusive prompt photon production. 
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polarization vector of a gluon is no longer purely transversal, but acquires an admixture of 
longitudinal and time-like components. Other important properties of the fc^-factorization 
formalism are the additional contribution to the cross sections due to the integration over 
the region above fi"^ and the broadening of the transverse momentum distributions due 
to extra transverse momentum of the colliding partons. 

In this approach, the treatment of the fcr-enhancement suggests a modification of the 
simple kT smearing picture described above: the transverse momentum kT of the incoming 
partons is generated in the course of non-coUinear parton evolution under control of the cor- 
responding evolution equation. First calculations of the inclusive prompt photon production 
at the Tevatron within the fc^— factorization formalism have been performed in [22-25]. The 
calculations [22-24] were based on the q + g* — > 7 + 5' and q + q J + g off-shell matrix 
element^. A reasonable agreement was found [24] between the theoretical predictions and 
the available D0 and CDF experimental data [1-6] in both the central and forward pseudo- 
rapidity 1]'^ regions. Perfect agreement was found also in the ratio of the cross sections 
calculated at ^/s = 630 GeV and y/s = 1800 GeV. However, an important component of the 
calculations [22-24] is the unintegrated quark distribution in a proton which at present is 
available in the framework of Kimber-Martin-Ryskin (KMR) [31] approach onljfl. In contrast 
with [22-24], the central part of our previous consideration [25] is the off-shell gluon-gluon 
fusion subprocess g* + g* — > 7 + qq. At the price of considering the 2 — > 3 rather than 
2 — s> 2 matrix elements, the problem of unknown unintegrated quark distributions has been 
reduced to the problem of gluon distributions. In this way, since the gluons are only respon- 
sible for the appearance of the sea but not valence quarks, the contribution from the valence 
quarks should be calculated separately. Having in mind that the valence quarks are only 
important at large x, where the traditional DGLAP evolution is accurate and reliable, this 
contribution has been calculated within the usual coUinear scheme based on 2 ^ 2 partonic 
subprocesses and on-shell parton densities. Thus, the way proposed in [25] enables us with 
making comparisons between the different parton evolution schemes and parametrizations 
of parton densities, in contrast with previous calculations [22, 24] where such a comparison 
was not possible. It is important that the predictions [25] based on the off-shell gluon-gluon 
fusion matrix element g* + g* ^ 'y + qq and the KMR gluon density agree with the previous 
results [24] based on the 2^2 subprocesses. This can be regarded as an additional proof 
of the consistency of the proposed method. 

In the present paper we will apply the formalism [25] described above to investigate 
the prompt photon and associated heavy (charm and beauty) quark production at high 
energies. The experimental data on the 7 + c and 7 + 6 cross sections as a function of photon 
transverse momentum have been reported recently [9] by the CDF collaboration. Also, 
there are available data [7, 8] on the associated prompt photon and muon production at the 
Tevatron, where the final state muon originates from the semileptonic decay of a charmed or 
beauty quark. Both these measurements are sensitive to the physics beyound the Standard 
Model (SM), for example the production of excited quarks or gauge-mediated supersymmetry 
breaking (GMSB) with neutralinos radiatively decaying to gravitinos [33]. Therefore, it is 
necessary to have a realistic estimation of the associated 7 + Q or 7 + /x production cross 
sections within QCD. An additional motivation for our investigations is the fact that these 

■^In the calculations [22], the usual on-shell matrix elements were embedded in precise off-shell kinematics. 
''^Unintegrated quark density was also considered recently in [32] . 
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processes provide a direct probe of the off-shell matrix elements g* + g* ^ + qq since there 
is no contribution from valence quarks. In order to investigate the underlying dynamics in 
more detail, we study the angular correlations between the transverse momenta of the prompt 
photon and the final muon. These quantities are sensitive to the production mechanism and, 
also, powerful tests for the non-coUinear evolution [34]. 

The outline of our paper is following. In Section 2 we recall shortly the basic formulas 
of the /cr-factorization approach with a brief review of calculation steps. In Section 3 we 
present the numerical results of our calculations and a discussion. Section 4 contains our 
conclusions. 



2 Theoretical framework 



2.1 Kinematics 

First, we recall in brief some technical details of our previous paper [25] needed below. 
We start from the kinematics (see Fig. 1). Let p^^^ and p^'^^ be the four- momenta of the 
incoming protons and p the four-momentum of the produced photon. The initial off-shell 
gluons have the four-momenta ki and k2 and the final quark Q and antiquark Q have the 
four- momenta pi and p2 and the mass uiq, respectively. In the pp center-of-mass frame we 
can write 

p(^^ = ^(1,0,0,1), p(^) = ^(1,0,0,-1), (1) 

where y/s is the total energy of the process under consideration and we neglect the masses 
of the incoming protons. The initial gluon four-momenta in the high energy limit can be 
written as 

ki = Xip(^) + kiT: k2 = X2P^^^ + k2T: (2) 

where kir and /c2r are the transverse four-momenta. It is important that kfj, — —k\rp ^ 
and \2^rj, — —k2T 7^ 0. From the conservation laws we obtain the following relations: 

kiT + = Pit + P2T + Pt, 

xiy/s = miTfi^' + m2Te^^ + Iprje^, (3) 

where y is the rapidity of produced photon, piT and p2T are the transverse four-momenta 
of final quark and antiquark, yi, 1/2, ^it and m2T are their center-of-mass rapidities and 
transverse masses, i.e. mfrp — mg -\- p|p- 



2.2 Cross section for associated 7 + Q production 

In general, according to the /cr-factorization theorem, the photon-quark associated pro- 
duction cross section can be written as a convolution 



J X2 27r 



a{p + p^-f + Q + X)= [ —A{xi,]4T,ii^)d\4j.^x 

J xi 27r 
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where a{g*+g* 1+QQ) is the partonic cross section, A{x, k|,, ij?) is the unintegrated gluon 
distribution in a proton and 0i and 02 are the azimuthal angles of the incoming gluons. The 
multiparticle phase space Iid?pi/2Ei5^'^^{J2,P'^~Yl,P°^^) is parametrized in terms of transverse 
momenta, rapidities and azimuthal angles: 



IT 



2Ei 2 " 2n 

Using the expressions (4) and (5) we obtain the master formula: 



dVirdyi—- (5) 



2 2^ A! l2 2Nn2 n2 > 2 2 > > > rf02 C^V^l #2 



(6) 



27r 27r 27r 27r 



where |A1p(5f* + f?* ^ 7 + QQ) is the off-mass shell matrix element squared and averaged 
over the initial gluon polarizations and colors, ipi and '?/'2 are the azimuthal angles of the 
final state quark and antiquark, respectively. Concerning the amplitude g* + g* ^ '~i + QQ, 
there are eight Feynman diagrams which describe this partonic subprocess at the leading 
order in and a (see Fig. 2). The analytic expression for the + 7 + QQ) 

has been derived in our previous paper [25]. We only mention here that, in accord with the 
/cT-factorization prescription [26, 27], the off-shell gluon spin density matrix has been taken 
in the form 

j:^^ihy''ik,) = ^- (7) 

In all other respects our calculations follow the standard Feynman rules. If we average the 
expression (6) over 0i and 02 and take the limit k^j. and h'^j' 0, then we recover the 
relevant expression in the usual collinear approximation. 

The multidimensional integration in (6) has been performed by the means of Monte Carlo 
technique, using the routine Vegas [35]. The full C++ code is available from the authors 
on requeslE]. This code is practically identical to that used in [25], with exception that now 
we apply it to calculate the cross section of prompt photon and heavy quark (or rather decay 
muon) associated production. 

2.3 Photon isolation and fragmentation contribution 

In order to reduce huge background from the secondary photons produced by the decays 
of 7r° and t] mesons, the isolation criterion is introduced in the experimental analyses. This 
criterion is the following. A photon is isolated if the amount of hadronic transverse energy 
E^'^'^ deposited inside a cone with aperture R centered around the photon direction in the 
pseudo-rapidity and azimuthal angle plane is smaller than some value E™^^: 

TTihad ^ T-imax 

< 5 /ON 

(^had_^)2^(^had_^^2<^2_ W 



^ lipatov@ theory, sinp . msu.ru 
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The CDF collaboration takes R = 0.4 and E^^^ = 1 GeV in the experiment [8, 9] and 
R = 0.7 and E^^^ = 2 GeV in the earlier experiment [7]. 

It is important that there is an additional mechanism of photon production not described 
above. It is the fragmentation of a partonic jet into a single photon carrying a large fraction z 
of the jet energy [8]. These processes are described in terms of quark-to-photon Dq^^{z, fi"^) 
and gluon-to-photon Dg^^{z, fj!^) fragmentation functions. However, the isolation condi- 
tion (8) not only reduces the background but also significantly reduces the fragmentation 
components. It was shown [36] that after applying the isolation cut the contribution from 
the fragmentation subprocesses is strongly suppressed (this contribution amounts to about 
10% of the visible cross section). Therefore in further analysis we will not consider the 
fragmentation component. 

3 Numerical results 

3.1 Theoretical uncertainties 

There are several parameters which determine the overall normalization factor of the cross 
section (6): the unintegrated gluon distribution in a proton ^(x, k^,/i^), the factorization 
and renormalization scales fip and fi^i and the heavy quark mass ttiq. 

Concerning the unintegrated gluon densities in a proton, we have tried here two different 
sets of them. These sets are widely discussed in the literature (see, for example, review [37] 
for more information). Here we only shortly discuss their characteristic properties. 

One set has been obtained [38] recently from the numerical solution of the CCFM equa- 
tion. Function ^(x,k|.,/i^) is determined by a convolution of the non-perturbative starting 
distribution Ao{x) and the CCFM evolution kernel denoted by A{x,'k'^, fi"^): 



In the perturbative evolution the gluon splitting function Pgg{z) including non-singular terms 
(as it was described in [39]) is applied. The input parameters in Aq{x) were fitted to 
reproduce the proton structure functions F2(x, Q^). An acceptable fit to the measured F2 
values was obtained [38] with x^/nrf/ = 1.83 using statistical and uncorrelated systematic 
uncertainties (compare to 'X^ /ndf ~ 1.5 in the coUinear approach at NLO). 

Another set (the so-called KMR distribution) is the one which was originally proposed 
in [31]. The KMR approach is a formalism to construct unintegrated gluon distribution from 
the known conventional parton (quark and gluon) densities. It accounts for the angular- 
ordering (which comes from the coherence effects in gluon emission) as well as the main 
part of the coUinear higher-order QCD corrections. The key observation here is that the 
dependence of the unintegrated parton distribution enters at the last step of the evolution 
and therefore single scale evolution equations can be used up to this stepcl. 

Significant theoretical uncertainties are connected with the choice of the factorization and 
renormalization scales. The first of them is related to the evolution of the gluon distributions, 

^In the numerical calculations we have used the standard GRV (LO) paramctrizations [40] of the coUinear 
quark and gluon densities. 




(20) 
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the other is responsible for the strong couphng constant As it is often done, we choose 

the renormahzation and factorization scales to be equal: fiR = fip = = ^\pt\- In order to 
investigate the scale dependence of our results we will vary the scale parameter ^ between 
1/2 and 2 about the default value ^ = 1. 

In the numerical calculations we set the charm and beauty quark masses to = 1.4 GeV 
and rub = 4.75 GeV. We have checked that the uncertainties which come from these quantities 
are negligible in comparison with the uncertainties connected with the unintegrated gluon 
distributions. For completeness, we use the LO formula for the strong coupling constant 
with Uf = 4 active quark flavors at Aqcd = 200 MeV (so that as(M|) = 0.1232). 
Note that we use the special choice Aqcd = 130 MeV in the case of CCFM gluon (as(M|) = 
0.1187), as it was originally proposed in [38]. 

3.2 Associated 7 + Q production at Tevatron 

We are now in a position to present our numerical results. Figs. 3 and 4 confront the 
7 + c and 7 + 6 production cross sections calculated as a function of the photon transverse 
energy with the preliminary experimental data [9] taken by the CDF collaboration at 
the Tevatron Run II (^/s = 1960 GeV). These data refer to the central kinematic region 
defined by \ri'^\ < 1. The solid and dotted histograms correspond to the results obtained 
with the CCFM and KMR unintegrated gluon densities, respectively. The upper and lower 
dashed histograms correspond to the scale variations in CCFM density as it was described 
above. We find that the CCFM-evolved unintegrated gluon density reproduces well the data 
within the theoretical and experimental uncertainties, and that the KMR density tends to 
underestimate the data in wide Et range. A similar effect was observed also in the case 
of inclusive prompt photon production [25]. The difference between the CCFM and KMR 
predictions is directly connected with the small-x behaviour of these gluon densities and 
demonstrates the importance of leading Inl/x terms. Of course, the scale uncertainties of 
our predictions are significant. The latter can be reduced by considering the ratio of the 
cross sections 7 + c to 7 + 6. This ratio is a subject of special interest: one could expect 
from the ratio of the quark charges that the 7 + c events should be 4 times more often than 
the 7 + 6 events. In addition to that, there must be extra suppression of the 7 + 6 events 
due to heavier b mass (in the g + g ^ 'J + QQ approach) or smaller beauty content in the 
proton sea (in the g + Q ^ J + Q approach). Our prediction for the ratio (7(7 + c)/(t(7 + b) 
is shown in Fig. 5 in comparison with the CDF data [9] . Both the CCFM and KMR gluon 
densities predict this ratio to be equal to 6 : 1 or 7 : 1 in a wide E^ range. This result is 
consistent with the measurement [9]. 

There are also available CDF data [7, 8] on the muons which originate from the semilep- 
tonic decays of charmed or beauty quarks. The experimental data [8] refer to the kinematic 
region pj. > 4 GeV, \r]'^\ < 0.9, \r]^^\ < 1.0 and y/s = 1800 GeV. In Fig. 6 we show the 7 + /i 
cross section as a function of photon transverse momentum pt- The contributions from both 
the 7 + c and 7+6 events have been taken into account. To produce muons from charmed and 
beauty quarks in our theoretical calculations, we first convert them into a D or B hadrons 
using the Peterson fragmentation function [41] and then simulate their semileptonic decay 
according to the standard electroweak theorjlJ. As usual, we set the fragmentation param- 

^Of course, the muon transverse momentum spectra are sensitive to the fragmentation functions. However, 
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eter e to Cc = 0.06 and = 0.006. The branching fractions /(c fi) and f{b fi) were 
set to /(c ^ ji) = 0.09 and f{h n) = 0.1078 [42]. One can see that in the case of 7 + 
production, our predictions with the CCFM gluon density shghtly overestimate the data but 
still agree with them within the uncertainties. The coUinear NLO QCD calculations [43] 
give similar description of the data. The results obtained with the KMR density lie below 
the measurements and are similar to those [8] obtained from the Pythia Monte Carlo [44]. 

An important point of our consideration is the investigation of the angular correlations 
between the prompt photon and heavy quark. It is well known that studying these cor- 
relations gives additional insight into the production dynamics and, in particular, into the 
effective contributions from higher-order QCD processes. For example, the lowest-order 
QCD production diagrams g + Q 1 + Q contain only the photon 7 and heavy quark 
Q in the final state. Therefore, the distribution over A0 = 0-y — 0q must be simply a 
delta function (5(A0 — tt) since the produced particles are back-to-back in the transverse 
plane and are balanced in pr due to momentum conservation. When higher-order QCD 
processes are considered, the presence of additional quarks and/or gluons in the final state 
allows the A0 distribution to be more spread and the heavy quark transverse momenta more 
asymmetric. In the /cT-factorization formalism, taking into account the non-vanishing initial 
gluon transverse momentum leads to the violation of back-to-back kinematics even at 
leading order. However, using the 2 — > 3 matrix elements instead the 2 — > 2 ones (as it 
was described above) makes the difference between the /cy-factorization predictions and the 
coUinear approximation of QCD (in aemC^'i approximation) not well pronounced. 

The associated 7+yU cross section as a function of the azimuthal angle A(6(7 — /i) has been 
measured in [7]. The data on the normalized differential cross section (1/a) da/dA(f){^ — ^i) 
have been presented. These data refer to the kinematic region YI < Et < 40 GeV, > 
4 GeV, \rf\ < 0.9, jr^^l < 1.0 and y/s — 1800 GeV. Our theoretical prediction compared to the 
data are shown in Fig. 7. We find here a number of the interesting points. First, the shapes 
of histograms predicted by the CCFM and KMR gluon densities are very different from 
each other. The CCFM gluon reproduces well the shape of the measured A0 distribution, 
although tends to shghtly overestimate the data at A0 ~ tt, while the KMR gluon density 
is unable to describe the data anywhere. The observed shape difference is in contrast with 
the transverse momentum spectra, where both the unintegrated gluon distributions under 
consideration demonstrate a (more or less) similar behaviour. This fact clearly indicates 
that the 7 + cross section as a function of A0 is very sensitive to the details of the non- 
coUinear evolution. A similar observation has been made earlier [34] in the case of 6-quark 
hadroproduction at the Tevatron. Thus, futher theoretical and experimental studying such 
quantities can give us the possibility to additional constrain the unintegrated gluon densities. 
However, we should mention that the behaviour of the A0 distribution at A0 ~ depends 
sensibly on the photon isolation criteria. In particular, it depends on the parameters R 
and determining the cone isolation (8). Our predictions for LHC energy are shown in 
Figs. 8 and 9. 

In conclusion, we would like to stress a number of important achievements shown by 
the /cj^-factorization approach. As a general feature, the model behaviour is found to be 
perfectly compatible with the available data on the heavy quark production as well as on 

this dependence is expected to be small as compared with the uncertainties coming from the unintegrated 
gluon densities in a proton. 
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the production of prompt photons and various quarkonium states at modern colhders [45]. 
It is important that the /c-r-factorization approach succeeds in describing the polarization 
phenomena observed in both pp and ep interactions (see, for example, [46] and references 
therein). The underlying physics is essentially related to the initial gluon off-shellness, which 
dominates the gluon polarization properties and has a considerable impact on the kinematics. 
So, we believe that the /cT-factorization formalism holds a possible key to understanding the 
production dynamics at high energies. Finally, once again we would like to point out the 
fundamental role of angular correlations which can serve as an important and crucial test 
discriminating the different non-coUinear evolution schemes. 

4 Conclusions 

We have tried a theoretical approach proposed in [25] to the associated production of 
prompt photons and heavy (charmed or beauty) quark in hadronic collisions at high en- 
ergies. Our approach is based on the /crfactorization scheme, which, unlike many early 
calculations [13, 16], provides solid theoretical grounds for adequately taking into account 
the effects of initial parton momentum. The central part of our consideration is the off-shell 
gluon-gluon fusion subprocess g* + g* ^ 7 + qq. At the price of considering the 2^3 
rather than 2 — 2 matrix elements, we have reduced the problem of unknown unintegrated 
quark distributions to the problem of gluon distributions. This way enables us with making 
comparisons between the different parton evolution schemes and parametrizations of parton 
densities, in contrast with previous calculations [22-24] where such a comparison was not 
possible (for the lack of unintegrated quark distributions except KMR). 

We have calculated the total and differential 7 + Q and 7 + yU cross sections (where muon 
originates from the semileptonic decay of the heavy quark Q) and have made comparisons to 
the recent CDF experimental data. In the numerical analysis we have used the unintegrated 
gluon densities obtained from the CCFM evolution equation and from the KMR prescription. 
It was demonstrated that the CCFM-evolved gluon density reproduces the Tevatron data 
very well, whereas the KMR gluon density is in disagreement with them. We especially point 
out the fundamental role of the angular correlations between the particles in the final state. 
These quantities can serve as a crucial test for the unintegrated gluon densities in a proton. 
Finally, we extrapolate the theoretical predictions to LHC energie^. 
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Figure 1: Kinematics of the g* + g* ^ ^ + QQ process. 
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Figure 2: Feynman diagrams which describe the partonic subprocess 5'* + 5'* — > 7 + QQ at 
the leading order in ag and a. 
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Figure 3: The cr{'~f + c) cross section as a function of photon transverse energy Et at \ri\ < 1.0 
and y/s = 1960 GeV. The sohd and dotted histograms correspond to the CCFM and KMR 
gluon densities, respectively, with the default scale /j, = Et- The upper and lower dashed 
histograms correspond to the scale variation in the CCFM distribution. The experimental 
data are from CDF [9]. 
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Figure 4: The (7(7 + 6) cross section as a function of photon transverse energy Et at \r]\ < 1.0 
and = 1960 GeV. Notations of histograms are the same as in Fig. 3. The experimental 
data are from CDF [9]. 
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Figure 5: The ratio of 7 + c to 7 + 6 cross sections as a function of photon transverse 
energy Et at |?7| < 1.0 and ^/s = 1960 GeV. The sohd and dotted histograms correspond 
to the CCFM and KMR gluon densities, respectively, with the default scale ji — Et- The 
experimental data are from CDF [9]. 
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Figure 6: The differential cross section da/dpT for associated prompt plioton and muon 
liadroproduction calculated at > 4 GeV, \r]'^\ < 0.9, \r]^\ < 1.0 and ^/s = 1800 GeV. 
Notations of histograms are the same as in Fig. 3. The experimental data are from CDF [8] . 
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Figure 7: The difference in azimuthal angle A0(7 — /i) between the photon and muon 
calculated at 17 < pt < 40 GeV, > 4 GeV, \r]^\ < 0.9, \r]^\ < 1.0 and = 1800 GeV. 
Notations of histograms are the same as in Fig. 5. The experimental data are from CDF [7]. 
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Figure 8: The differential cross section da/dpr for associated prompt photon and muon 
hadropro duct ion calculated at > 4 GeV, Irj"'] < 2.5, \7]^\ < 2.5 and a/s = 14 TeV. 
Notations of histograms are the same as in Fig. 5. The isolation criterion as in [8, 9] was 
applied. 
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Figure 9: The difference in azimuthal angle A0(7 — //) between the photon and muon 
calculated at 10 < < 100 GeV, > 4 GeV, \ri^\ < 2.5, Ir/^^l < 2.5 and ^ = U TeV. 
Notations of histograms are the same as in Fig. 5. The isolation criterion as in [8, 9] was 
applied. 
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